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Abstract 

We critically analyze recent attempts to include final state interaction (FSI) ef- 
fects in the calculation of the CP violating ratio e'/e. In particular the approach 
of using the Muskhelishvili-Omnes equation for the decay amplitudes is examined. 
We find that the values of the dispersive correction factors are very sensitive to 
initial conditions which are, at present, out of control, and demonstrate that the 
claimed large enhancement of e'/e through FSI is questionable. We propose, in- 
stead, a different implementation of this approach which may be useful to lattice 
determinations of weak amplitudes. 



1 Introduction 



One of the crucial issues for our understanding of CP violation is the question whether 
the size of the observed direct CP violation in Kl — > tttt decays, expressed in terms of 
the ratio e'/e, can be described within the Standard Model. Experimentally the grand 
average, including the results from NA31 0, E731 §, KTeV § and NA48 f§ reads 

Re(e'/e) = (21.2 ± 4.6) x 10~ 4 . (1) 

There are different opinions whether this result can be accommodated within the 
Standard Model. In [(J the central value of the predicted e'/e is typically by a factor of 
2-3 below the data, with an estimated theoretical uncertainty of the order of 60 -r- 100 %. 
The result in Eq. (fj), then, can only be accommodated if all relevant parameters are 
chosen simultaneously close to their extreme values. Higher values of e'/e, compatible 
with the data within one standard deviation have been found, instead, in [0, §. Recent 
reviews can be found in |J. 

More recently, following the previous work of Truong [Til, it has been pointed out by 



Pallante and Pich |TlJ that the inclusion of final state interactions substantially affects the 
estimates of e'/e. The size of the effect is parametrized in terms of dispersive correction 
factors IZi multiplying the corresponding isospin amplitudes Aj. With TZq « 1.4 and 
7^-2 ~ 0.9, as found in their numerical evaluation, a substantial enhancement of e'/e was 



obtained. Indeed, by including these factors, the authors of Ref. j[lj find that the central 
value e'/e = 7.0 x 10 -4 of |J gets increased to ~ 15 x 10~ 4 which is much closer to the 
experimental average ([I]). Similar results, although with a different dispersive approach, 
have been found in JTJJ . 

The interesting paper of Pallante and Pich stimulated us to give a closer look to the 
whole procedure followed to obtain the correction factors TZj. The phenomenological 
application of the Muskhelishvili-Omnes equation to the K — > -kit case requires some 
arguable assumptions (see Section H). However, even within those assumptions, we found 



a major objection which unfortunately makes the calculation of Ref. [IT] questionable 



On the basis of symmetry arguments only, one can show that there is a point, out of 
the physical cut, where the amplitudes induced by left-handed operators are expected to 
vanish. As demonstrated below, the authors of Ref. [|TTJ] made an assumption on the value 



of the derivative of the amplitude at this point which is not justified. All their results 
rely on this assumption. We will show that, by replacing the latter condition with others, 
equally acceptable on physical grounds, it is possible to obtain very different values of 
the correction factors 1Zi. Thus we conclude that the results found in |T| are subject to 
substantial uncertainties. 
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In spite of this, we found some interesting application of the dispersive analysis to 
non-perturbative calculations of weak amplitudes. For example, lattice calculations could 
provide a condition which replaces the arbitrary assumption on the derivative of the 
amplitude and fix unambiguously the factors IZi (under the hypotheses inherent to the 
use of the Muskhelishvili-Omnes equation in K — > tttt decays). 

Our paper is organized as follows. In section 2 we present general formulae for the 
Muskhelishvili-Omnes solution of dispersion relations applied to K —>■ tttt amplitudes, 
focusing on the assumptions needed to obtain the final result. Using these formulae we 
demonstrate in Section 3 that the evaluation of the dispersive factors TZi along the lines 



proposed in |T]J is ambiguous. In particular we stress that there is no justification to use 
the numerical values of TZj found in [|11[] in conjunction with present lattice and large- Nq 
calculations of hadronic matrix elements. Similar comments apply to |12]. In Section 
4 we briefly illustrate how modifications of the technique developed in [10, [ll|] could in 
principle improve the accuracy of the predictions for e'/e in conjunction with future lattice 
calculations. We conclude in Section 5. 



2 General formulae 

In this section we recall the basic ingredients used to derive the solution of the Muskheli- 
shvili-Omnes equation [T3|] in the case of K — > tttt amplitudes. The starting point is the 
iV-subtracted dispersion relation. We denote by A(s) the decay amplitude of a kaon of 
mass ttlk = y/s into two pions of mass m n . The invariant mass of the two pions is taken as 
(pi +P2) 2 — s, corresponding to the insertion of a weak operator carrying zero momentum. 
The following discussion applies as well to the case where A(s) is the K — > tttt matrix 
element of a local operator renormalized at a fixed scale \x. 

Following Ref. [TTJ we assume that A(s) is analytic in the cut s plane, with the cut 
going from 4m^. to 00 []. Further assuming that s~ N A(s) — ► for s — ► 00, we can write 

. , « „ , , (s-s ) N r°° , ImA(z) , , 

A(S) = P„_ l( .,. o) + ^T~jim> ~(z - s^ijls - Q ' 

where Pjv-i(s, so) is a — 1 degree s polynomial, the coefficients of which are to be fixed 
by iV independent conditions on the amplitude and so denotes the subtraction point (in 
general one could choose N different subtraction points), which has to lie outside the 
physical cut. Writing the imaginary part of the amplitude as 

lmA{s) = A{s)e- iS(s) sin 5{s), (3) 



1 These analiticity properties of A(s) are certainly not fully correct. However one can argue that 
Eq. (||) still holds to a reasonable accuracy in a region close to the physical point (see Section 4) . 
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the dispersion relation (|2]) becomes a Muskhelishvili-Omnes equation |IJJ. Here 5(s) is 
the strong phase of the amplitude, which in the elastic region is equal to the phase shift 
of the tttc scattering appropriate to a given isospin channel Q. The general solution of this 
equation is given by 

A(s) = Q N - 1 (s)0( 8 ), (4) 

where 

0(s) = exp (- r dz . 5{Z) ) (5) 

\7T J Ami [Z — S — ie) J 

and Qat_i(s) is a iV — 1 degree polynomial. Equation (|) shows explicitly that N indepen- 
dent conditions on A(s) and the knowledge of the phase on the whole cut are sufficient 
to fully determine the amplitude. Note that there is no longer any reference to the sub- 
traction point So of the original dispersion relation. 

We point out that, if we were able to obtain a number of conditions on A(s) larger 
than N, we could reduce the sensitivity of the solution to the knowledge of 5(s) at large 
s. Indeed, we can always rewrite the solution (f|) as 



A(s) = Q N -i(s) exp (^J2 c i s ^j 



7T 



where the product of the first two terms is completely determined by iV + M independent 
conditions. If possible, this could be useful since, in practice, S(s) can be extracted from 
the data only up to the inelastic threshold. 

To make contact with the work of Pallante and Pich [|Il|] , we now discuss two examples 
on how A(s) is determined from a specific set of conditions in the case N = 2. In this 
case, corresponding to an amplitude going at most linearly in s for large s, two conditions 
on the amplitude are required: 

1. In order to fix the coefficients of Qi(s), a possible set of conditions is given by the 
knowledge of the amplitudes at two different points si and S2, namely 

A(s 1 ) = A 1 , A(s 2 ) = A 2 . (7) 

In this case, one finds 

A(s)= — w^moM — 0(s) - (8) 



2 Note that even in the elastic region S(s) can be identified with the phase shift measured from tttt 
scattering only if the dependence of the latter on the kaon mass is negligible. We thank G. Colangelo for 
pointing this out to us. 
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We stress that, although the r.h.s of this equation is written in terms of si^, the 
physical amplitude is obviously independent of the choice of these points, provided 
that Ai ; 2 are given by Eq. (0). 

2. A different set of conditions can be given by 

A(si)=0, A\s x ) = C, (9) 

namely by knowing the point s\ where the amplitude vanishes, or assumes a specific 
value, and its first derivative in that point. With the conditions in Eq. (§) the 
solution is particularly simple and takes the form 

A(s) =C(s-s 1 )0(s 1 )- 1 0(s). (10) 



3 Why we are not able to compute universal disper- 
sive factors 



We now critically examine the procedure followed by Pallante and Pich |TT| to compute 



the so-called dispersive factors TZi for the 1 = and 1 = 2 amplitudes. To this end, we 
start by summarizing the main steps of the procedure adopted in Ref. JTT . 



On the basis of the lowest-order chiral Lagrangian, Pallante and Pich assume that the 
solution of the Muskhelishvili-Omnes equation is given by 

A I (s)=C I (s-ml)n I (s), (11) 

where 

Here Oj is given by Eq. fl5|) with 5 replaced by 5i, the experimental elastic irir(I) phase 
shift, integrated up to a certain cutoff. This solution corresponds to the case in Eq. fllOD 
with si = m%, i.e. the authors of Ref. [Til] impose Ai(m%) = and implicitly assume the 



knowledge of A'^rnD = In order to fulfill the latter condition, they assume without 
justification that the matrix element at the physical point, 

M I = (nn(I)\H cS \K)\ s=m2K , (13) 

provided by non-perturbative methods such as lattice QCD or large iV^ estimates, corre- 
sponds to C/(s — m%)\ s=m 2 . Their solution can thus be written as 

Ms) = M T { "- m l ] e^TZ^s) (14) 

\ m K _ m TT> 



The need of these two assumptions is common also to Ref. jL0[, as pointed out in 114 
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and the dispersive factors TZi are simply given by 

Hi = ni{m 2 K ) 







Oi(m\) 


Mi 




Oi(ml) 



(15) 



We have the following comments on this approach: 



One of the two conditions used to determine the solution is Ai(m 2 ) = 0. This is 
justified in Ref. |TT]] on the basis of the lowest-order chiral realization of (8^, 1#) 
and (27 l, 1r) operators. The zero at s = m 2 of the amplitudes induced by these 
operators actually holds beyond the lowest order |15|, [16|] . Indeed, it is a consequence 
of the vanishing of these amplitudes in the SU(3) limit, as follows by an old result 
by Cabibbo and Gell-Mann [[17]]. In modern language, we can rephrase the result of 
Ref. [I7j by saying that the on-shell K — > 2n matrix element of any local (8l, li?) 
operator invariant under CPS symmetry vanishes in the absence of SU(3) 
breaking. Therefore we think that this condition is well justified in the case of 



[8l, 1r) operators and, particularly, in the case of Q e considered in HTT | . 



Our main criticism is related to the determination of the second condition on Aj(s), 
namely the one on A'^m 2 ). Chiral perturbation theory alone cannot fix the value 
of this derivative, which must be known from some non-perturbative information. 
In their procedure, Pallante and Pich implicitly assume that the matrix elements 
determined by non-perturbative methods such as lattice or large Nc, fix the value 
of A'j(m 2 ) according to the condition 



A'i(s 



Ci 



Mi 



m 



K 



mi 



(16) 



This assumption, which may look reasonable on the basis of the lowest-order chi- 
ral Lagrangian, actually involves an ambiguity of the same order of the dispersive 
correction itself. To illustrate this point, we make a different choice of the initial con- 
ditions (while keeping the same Mi): we use Mi to fix the value of the amplitude 
at threshold, via the relation 



Aj(s = Ami 



Mi 



3mz 



m\ — m\ 



(17) 

provides a good approx- 

,2 



In other words, we assume that Mi(s — m^)/{m 2 K — m 2 , 

imation to the real amplitude near s = Am 2 instead of the point s = m~ as was 
implicitly employed in [11]]. One may argue that this is also a reasonable choice 
because strong-interaction phases vanish at threshold, but of course this condition 
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is as arbitrary as the one adopted in [IT]. Using this condition and A(s 
from Eq. ([8]) we find 



ml 



Axis) = Mr 



(s-mDOjjs) 
m 2 K — rn%.)Oi(4m%) 



(18) 



With the same parameterization of the phase Sq(s) as in Ref. |TTJ and our choice 
(|I8"D, the dispersive factor for (8^, 1r) operators is 



threshold 




K) 



O (4ml) 



1.1 



instead of 



K 



pp 



O (m 2 



K) 



1.4. 



(19) 



(20) 



O (ml) 

With TZ = 1.1 (and TZ 2 = 1.0) the central value e'/e = 7.0 x 10~ 4 of |S| gets increased 
to ~ 9 x 10~ 4 , instead of ~ 15 x 10~ 4 found in |lTJ. While the first enhancement is 
well within the theoretical uncertainties quoted by the various analyses [§]-[§, the 
second one would have a considerable impact on the predicted value of e'/e. 

Other choices of the initial conditions, equally acceptable under similar assumptions, 
would lead to still different results for TZ . This exercise illustrates that, unless one 
knows the value of s at which a given non-perturbative method provides the correct 
value of the amplitude (and/or its derivative), it is impossible to unambiguously 
compute the appropriate correction factor due to final state interactions. 

In the case of AI = 3/2 matrix elements of (8l,8r) operators, such as (Q$ ) 



considered in Ref. the amplitude does not vanish at s = m^. Still the authors 
of [|TTJ claim that the correction factor is 1Z 2 = |0 2 (m|-)/02(m^)| ~ 0.9. We 
can reproduce this result assuming that the dispersion relation requires only one 
subtraction and that the non-perturbative calculations give the correct amplitude 
This condition, however, appears even more arbitrary than the one 



ml 



at s 

used to determine A'i^m 2 ) in the case of (8/,, 1r) operators. Indeed in this case the 
point s = m\ does not play any special role. It is therefore not surprising that in a 



different but consistent framework, as the one discussed in |L9], ^(J , non-perturbative 
corrections to (Q^ 2 ^) are found to enhance, rather than suppress, lattice and large- 
N c estimates. 

Let us then summarize the main point we want to make. While the full amplitude 
clearly does not depend on the choice of the value of s at which the conditions on the 
amplitude are placed, the dispersive correction factors are obviously dependent on this 
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choice as we explicitly demonstrated above. As presently we are not in a position to state 
precisely which value of s the exisiting non-perturbative approaches, like lattice, 1/N, 
etc., correspond to, it is not possible to determine uniquely the corresponding factors IZi 
to be incorporated in each of those methods. 

Before concluding this section, we comment on other attempts present in literature 



to include FSI effects in e'/e. The approach of Ref. |T2fl is somehow similar to the one 



of Ref. lyj. In this paper, however, the integral over the cut is divided in two parts: a 
low-energy region from threshold up to /i 2 and a high-energy region from /i 2 up to oo. By 
identifying the scale /x with the renormalization scale of the inserted operator, the high- 
energy region is claimed to be estimated in perturbation theory. The results thus obtained 
are in qualitative agreement with those of Ref. [Tlj] and show a rough independence from 
the unphysical scale \i. We stress, however, that the matching procedure adopted in this 
paper cannot be justified theoretically. Since we do not see a way of deriving consistently 
the formulae of Ref. [EJ from first principles, we are skeptical about the approach followed 
in this paper. 

Other recipes to include a posteriori FSI effects in non-perturbative calculations of 
weak amplitudes were presented in the recent literature J7|, [Hp. For these we have ob- 
jections similar to those made to [[11]]: it is not clear to us how to relate the theoretical 
values of the amplitudes, computed with some non-perturbative method, to the physical 
ones which include FSI. In cases where the results of non-perturbative calculations are 
unable to reproduce the experimental phases, it is not possible to extract unambiguously 
the absolute value of the amplitudes. Since corrections may be applied only if this issue is 
under control, we think these procedures to incorporate empirically correcting factors to 
account for FSI do not help to improve the accuracy of the calculation. These attempts 
provide at best estimates of the theoretical uncertainties, as for instance discussed in M . 



4 Further developments 

In spite of the criticism discussed before, we believe the Truong-Pallante-Pich approach 
is interesting and deserves further investigation. For this reason, in this section we briefly 
discuss possible applications of this technique to improve the accuracy of theoretical pre- 
dictions for K — > TT7C matrix elements. We first want to state the assumptions needed to 
make this approach useful: 

1. In general, in order to fix the number N of subtractions in the dispersion relation, 
hence the degree of the polynomial Qn-i, one has to make an ansatz on the be- 
haviour of the amplitude as s — > oo. In the present case, however, Eq. flj|) should 



7 



be understood as an approximate equation, valid in a limited s region close to the 
physical point. Indeed other sources of non-analiticity besides the cut starting at 
4m^ can appear in A(s) due to the identification of s with kaon mass, that is a 
parameter of the underlying theory. Still the solution of Eq. (^|) can give a good 
estimate of the true amplitude if these further sources of non-analiticity are well rep- 
resented by a polynomial expansion in the considered region. This is for instance 
the case in an effective theory where the pions are the only dynamical fields, or in 
the quenched approximation of QCD. In this context the number of subtractions is 
no longer related to the behaviour of the amplitude at large s, rather to the accu- 
racy one needs to achieve. Large values of N are obviously preferable, although in 
practice they are hardly attainable. Therefore, following Refs. [pi [IT], we assume 
that N = 2 provides a reasonable approximation. 

2. Above the inelastic threshold, 5(s) is no longer the tttt phase shift and cannot be 
directly extracted from the experimental data. In addition, even in the elastic 
region, the phase of A(s) can be identified with the experimental tctc phase shift 
only assuming that the dependence of the nir scattering amplitude on the kaon 
mass is negligible. Were we able to compute the weak amplitude for a large number 
of points, we could reduce the sensitivity of the solution to the value of 8(s) at 
large s, as shown in Eq. (||). As already stressed, this possibility is presently rather 
remote. In practice, the integral entering O(s) has to be restricted to the elastic 
region, unless some further assumptions on the amplitudes of the other channels are 
made. These issues entail some uncertainties which have been partially discussed in 



Ref. Ill 



Under these assumptions, we now envisage a possible implementation of this approach 
in lattice simulations. It has been shown that, even in the presence of chiral symmetry 
breaking and of the Maiani- Testa theorem lattice calculations can obtain the physical 
matrix elements at threshold, corresponding to s — Am^ [23]. This information, once 
available, combined with the constraint A(s = m^) = and with the ansatz that the 
amplitude grows linearly in s, is sufficient to determine the physical result for the real 
kaon mass (see Eq. (§)) in the case of (8l, 1r) operators. 

Concerning AI = 3/2 matrix elements of (8l,8r) operators, one may argue that one 
condition is not enough to fully determine the amplitude. We stress, however, that in this 
case the extrapolation between the tttt threshold and the physical kaon mass is expected 
to be smoother, resulting in a correction factor which can be approximated by unity with 
high accuracy. 
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The evaluation of the amplitude at s > 4m^ is particularly difficult with lattice simu- 
lations since these are performed in the Euclidean space-time. Proposals, however, exist 
for overcoming this problem f24] , |25| . We finally stress that the region s < 4m^ is not 
accessible to lattice calculations, as well as to any other approach, if the particles ought 
to be on-shell and the inserted operator carries zero momentum. 

5 Conclusions 

We have critically discussed the issue of final state interactions effects in the evaluation 
of the ratio e'/e. Whereas we cannot exclude possible substantial enhancements of e'/e 
through FSI over lattice and large-iVc estimates, our analysis demonstrates that the 
present calculation of the dispersive factors IZi is subject to considerable uncertainties, 
which prevent to draw any definite conclusion on the importance of these effects. 

We summarize the main point of this paper. While the full amplitude clearly does not 
depend on the choice of the value of s at which the conditions on the amplitude are placed, 
the dispersive correction factors are obviously dependent on this choice. As presently we 
are not in a position to state precisely which value of s the exisiting non-perturbative 
calculations correspond to, it is not possible to determine uniquely the corresponding 
factors TZi to be incorporated in each of those methods. 

On the other hand, we have shown that the Muskhelishvili-Omnes approach could 
be useful to reduce, or at least to better estimate, the uncertainties on the physical 
amplitudes, once lattice results on the matrix elements at threshold will be available. 
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